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In this paper we propose a general framework for deriving 
the effective interactions for many-body systems. We show 
how the selfconsistent interaction can be constructed on the 
quantum level for both local and nonlocal potentials. In par- 
ticular the relation between the selfconsistent and symmetry 
restoring interaction is derived. The difference between the 
two is related to the coupling of the perturbing field to the 
kinetic energy operator. 

PACS numbers: 21.60.-n, 21.60.Jz, 21.30.-x 

For the study of collective motion of the nuclear or 
other many-body systems it is essential to properly de- 
fine the residual interactions. This problem appears al- 
ways when the ground state has been obtained by means 
of some variational principle without the explicit use of 
many-body forces. Based on the requirement of selfcon- 
sistency, Bohr and Mottelson derived the effective in- 
teractions for spherical nuclei Q. For a deformed sys- 
tem, like the Nilsson potential, Marshalek has shown 
that the volume conservation constraint is equivalent to 
a Hartree approximation for many-body forces The 
effective interaction can be formally constructed as func- 
tional derivatives of the total energy with respect to the 
density. For the harmonic oscillator potential it can be 
done analytically H,^. Based on this method and the 
Thomas-Fermi approximation it has been shown how the 
subsequent terms of the effective interaction can be ob- 
tained |,|]. 

However, all these methods unnecessarily rely on semi- 
classical approximations. Moreover, the comparison be- 
tween the selfconsistent strengths of the interaction and 
the so-called microscopic ones based on the linear re- 
sponse theory can result in different values The na- 
ture of this difference, albeit attributed to shell effects, 
seemed not to be clear, but can be elucidated on the 
quantum level. Last not least, there is a powerful sym- 
metry restoration method introduced originally by Py- 
atov and collaborators Q and developed in a series of 
papers |^ , which appears to be related to the concept of 
selfconsistent forces. To our knowledge this link has not 
been established yet. In the present paper we show how 
one can formulate the problem on the quantum level. We 
will not address any specific choice of coordinates which 



could enable us to decouple collective modes from spuri- 
ous components. This can be done at any stage of the 
presented method, although it is not crucial for the merit 
of this article. 

Let us assume that our phenomenological one body 
hamiltonian has the form: 



H = T + V, 



(1) 



with the kinetic energy T — and the potential 

2m 

V{r,r'). The nonlocality of the potential V can appear 
e.g. as a consequence of a spin-orbit couphng. The so- 
lution of the hamiltonian (|l|) defines a one body density 

N 

p(r, r') = ^0fe(r)0fc(r'), where (j)k{r), k = 1,...,N de- 

fc=i 

note the N lowest eigenfunctions of H, and the density 

N 

operator is defined via p — ^~^|fc)(fc|. To lowest order 

the hamiltonian describing the selfconsistent interaction 
is assumed to have the separable form: 



1 -2 
Hint = -j^xQ , 



(2) 
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where Q is a hermitian one-body operator and x denotes 
a coupling constant. The reason that the particular form 
of the interaction (||) is used, instead of the more general: 

Q' Q' , with Q' 7^ Q', is associated with the fact that 
the exchange potential generated by the separable inter- 
action can be neglected. Hence the Hartree fields of the 
interaction Q' Q' and (H) are equivalent if Q is a her- 
mitian or antihermitian part of Q' . Also in the RPA ap- 
proach the interaction Q' Q' can always be transformed 
to the form (H) in the quasiboson approximation. 

The following conditions should be imposed: (i) 
Tr{Qp) — {Q)o = 0, (ii) the interaction Hint should 
induce the same change in the density p as in the poten- 
tial V . Clearly, condition (i) provides a restriction on Q, 
whereas (ii) determines the coupling constant. Condition 
(i) is rather difficult to fulfill for an arbitrary potential, 
except when V is spherical. For the deformed harmonic 
oscillator, Sakamoto and Kishimoto derived the double- 
stretched multipole interaction that always fulfills crite- 
rion (i) [Q. This interaction has the further advantage 
to exactly decouple the spurious K = 1 modes from the 
shape multipole oscillations. Recently Kubo et al. 0| 
extended these methods to momentum dependent po- 
tentials. For realistic potentials like the Woods-Saxon 
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potential, condition (i) can be fullfilled approximately in 
a similar manner [ p^ . 

To analize the role of the condition (ii), assume that 
the collective mode is produced by the one body field A. 
Hence, the operator U = ex-p{—iAd) defines the transfor- 
mation of the ground state (iV— body) wave function: 



(3) 



Since we limit ourselves to small perturbations we will 
consider only the first term of the expansion. In order 
to preserve the normalization of the wave function we 
require A to be hermitian {9 is real) . Instead of changing 
the wave function it is easier to consider the change of 
the operators V and p: 



U+VU KiV + i9[A,V] 
U+ pU K p + i0[A, p]. 



(4) 
(5) 



The condition (ii) requires that is proportional to 

{[in-- 

= ix{[A, V]) - *xTr([i, V]U+pU), (6) 
which leads to the equation: 



= tx{[A, V])o~xO{[[A,V],A])o. 



(7) 



Since in general the first term need not to vanish, an 
arbitrary field A does not generate a selfconsistent mode 
unless one will make the replacement: [A, V] [A, V] — 
{[A,V])o. Only when: ([i,f])o = 0, eq. (0) defines the 
selfconsistent interaction strength x^'^'''^ : 



X 



self 



'{[[A,V],A])o. 



(8) 



Thus the selfconsistency requirement specifies the inter- 
action: 

H,nt = lx'''^Q' = -lx''''^[A,V?- (9) 

Note, that if one chooses A as an operator associated 
with a particular symmetry like e.g. angular momentum 
component Ji, or total momentum P, then the above 
prescription is equivalent to the Pyatov formula for the 
symmetry restoring interaction ||^ . Since these operators 
commute with the kinetic energy T, one can write in 
the form: 



Hi, 



1 



with the symmetry restoring coupling strength: 
1 



X' 



= -{[[A,H],A])o. 



(10) 



(11) 



In order to consider excitations that cause a change of 
the density, we need to find an explicit form of the oper- 
ator A. Starting from the single-particle wave functions 
we note that they will be transformed according to: 



Mr) 



1 



VI - V • u 



fc(r + u(r)). 



(12) 



where u is in general a function of r. Since the ground 
state is an independent particle state, eq. (12) generates 
a transformation of the many body wave function ip. The 
Pauli principle is not violated by this transformation be- 
cause all the wave functions are transformed in the same 



manner 1 1 1 1 . The transformation (y_2|) is generated by an 
operator A of the form: 



1 



and 



A = 7:(u(r) • p + p • u(r)). 



0fc(r) + 0u.V0fe(r) + -^ 



(13) 



fc(r)V-u. (14) 



The field A contains e.g. all the 2'^-pole distortions. In- 



deed, by choosing u(r) = VQa^, where Qa^ 



1 



(Qam- 



) or I 



V2 



(Qam-Q 



V2' 

the wave function ip will 



be affected by the multipole field. 

Following the above prescription one can easily (at 
least formally) construct the appropriate selfconsistent 
interaction: 



Hint 

1 



^X'^'nu-p 



X 



self 



[up 



p u, \>],u p 



(15) 



P • u])c 



This interaction is not anymore limited to the volume 
conserving modes, V • u = 0. Other excitations like the 
breathing mode can also be considered. In the case of a 
multipole field and a potential which does not depend on 
the momentum the relations (^5|) reduces to the expres- 
sion: 



Hint 

1 



1 



X 



self 



x^'^'-^(VQAM-vy)2, 

-(V(VQam-VF)-VQam>o. 



(16) 
(17) 

because 



This result differs slightly from that in Ref. 
we consider only a single excitation mode. 

In the case of nonlocal potentials expression ( p^ is 
not valid since [VQau, ^ and one has to use the 
general expression (lla). Note, that the transformation 
( p^ ) induces a deformation of the Fermi sphere resulting 
in a change of the current distribution: 



jo + (u • V)jo + (jo • V)u + joV • u, (18) 
^V(/)fe — (V(/)^)(/)fc), and jo denotes the 



1 ^ 

where i = — > 
equilibrium distribution. Hence, an additional correction 
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arises for nonlocal potentials. This effect was not con- 
sidered in previous approaches since it vanishes in the 
hydrodynamical limit. It should be mentioned that since 
the nonlocal potentials violate the Galilean invariance its 
presence lead to the appearance of an additional interac- 
tion discussed in Ref. p^ . 

Let us now examine in more detail the concept of self- 
consistent and symmetry restoring interactions. If the 
hamiltonian of eq. (1) violates a given symmetry one may 
restore it in the RPA approximation by adding an inter- 
action of type (10). The interaction strength x*"*^^* given 
by the expression (11) assures that there are no restoring 
forces towards the perturbation generated by the field 
A 1^. For the harmonic oscillator model the equality 
between the latter and the selfconsistent strength holds 
only for a few particular cases whereas usually it is dif- 
ferent by the value of a discrete jump, attributed to the 
influence of the shell effects . In the following we ar- 
gue instead, that these jumps are related to two different 
levels of symmetry restoration. The selfconsistent inter- 
action arises from the requirement of the invariance of the 
potential towards the perturbation generated by field A, 
whereas the symmetry restoring forces assure the invari- 
ance of the total hamiltonian. 

As discussed above, the operator A does not commute 
in general with the kinetic energy operator. Hence, one 
cannot replace the potential V by the total hamiltonian 
H in formula (||). This is a consequence of the fact that 
in general the perturbed hamiltonian and density do not 
commute. Namely: 



LUiTii = const., 



(19) 



where H is a one body field generated by the density p 
through the Ilartree(-Fock) procedure i.e. [H{p),p] — 0. 

This can easily be verified in the case of the harmonic 
oscillator model. Consider the three-dimensional spheri- 
cal harmonic oscillator hamiltonian: 



H 



P 

2m 



(20) 



Assume that the operator A is choosen in the form: 

i = VQ2o-p. (21) 
Then the perturbed hamiltonian has the form: 

H' = H + 2muj'^eQ20 (r) . (22) 
If we define a = 49 then: 



iJ' = 1^ + yc^'((l - a)xl + (1 - a)xl + (1 + 2a)xl) 



-A2 



2™ 2 ^ 



2 2 
UJ, X., . 



(23) 



i=l,2,3 



For the harmonic oscillator the selfconsistency require- 
ment between the potential and the density takes the 
simple form: 



(24) 



where — mujiTr{xfp). However, if one calculates the 
eigenfunction of the perturbed hamiltonian then relation 
(p3) is violated, since 



where 



Sloi — 5lo2 — — —auj + O(a^) 
Slj^ — auo + 0{a^). 



(25) 



(26) 



We assumed that there are no configuration changes 
caused by the perturbation which obviously is the case 
for the small perturbations considered in this article. 

To preserve the selfconsistency in the perturbed sys- 
tem one has to transform the wave function according 
to eq. (3), where the operator A is defined by eq. (21). 
It defines the new single particle wave functions 4>'f. (r) — 
4>k{Y'), where 



r' = r+-VQ20. 
Using such wave functions one can easily find that: 



ujiT,^ = uj^^{nk + -) ^ const.. 

k=l 



(27) 



(28) 



and the relation ( p^ ) is preserved. 

Hence we see that two methods of obtaining the per- 
turbed wave function do not coincide. Only one of them 
generates the selfconsistent density of the perturbed sys- 
tem. This difference is responsible for the difference be- 
tween the selfconsistent and symmetry restoring coupling 
constants. To make the discussion more general we use 
now the response function formalism to show when these 
two concepts lead to the same results. 

Consider the perturbation of the hamiltonian defined 
as: 

SH = i[A, H] = i[A, f] + i[A, V]=5f + 5V. (29) 
The total hamiltonian has the form 



H' = H + dH0{t), 



(30) 



where the perturbation in general can be time-dependent. 
We define the quantity: 



AE{t) = {H') - (i/')o, 



(31) 



which represent the difference between the expectation 
value of the hamiltonian in the perturbed and initial 
state. The Fourier transform of this quantity is related 
to the Fourier transform of the function 9{t): 



AE{u)=RsH,5H{L0)e[L0). 



(32) 
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The response function Rsh.sh{^) has the form: 



This impHes that: 



Rsh,sh{i^) = - 

fe>o 



\{k\5H\Q)\HEk ^ Eg) 
{Ek - Eor - to^ ' 



(33) 



ST = 0. 



(41) 



where Eq and Ek denote the energies of the ground state 
|0) and excited states |fc), respectively. 

On the other hand the evaluation of the left-hand side 
of Eq. ^ atLu = gives: 



AEiu;^0)^-([[A,H],A])o. 
Hence in the case when [T, A] = we get: 
1 



X' 



self 



Rsh,sh{ijJ = 0) = R&v.sv{i~o = 0), 



(34) 



(35) 



where x"^^^ has been defined according to (^) . This result 
is exactly the prescription for determining the symmetry 
restoring coupling strength. One can see that it coin- 
cides with the selfconsistent strength when the operator 
A commutes with the kinetic energy operator. 

In the case when [T, A] 7^ one can rewrite the equa- 
tion (B2) in the form: 



Consequently for the dipole mode the microscopic and 
selfconsistent strength coincide: Xijj''* = xT^l^ ■ 

To conclude, in the present paper we have shown that 
the selfconsistent interactions can be derived at the quan- 
tum level for general potentials. The Pyatov prescrip- 
tion for the symmetry restoring forces can be included in 
this general framework. The change of the momentum 
distribution will produce an additional correction to the 
selfconsistent interaction if momentum dependent poten- 
tials are considered. The origin of the difference between 
the symmetry restoring and selfconsistent coupling con- 
stants relates to two different levels of symmetry restora- 
tion. Whenever the perturbing field commutes with the 
kinetic energy operator the two concepts are equal. 

PM acknowledge the financial support from the 
Swedish Institute and the Goran Gustafsson foundation. 



1 



.self 



R5V,5v{i^ = 0) + Rsv,5TiuJ = 0), 



(36) 



R&v.st{^ = ^) = -y^\ ? ? + c-c- , (37) 



fe>C 



where we have used the relation: 



{[[A,f],A])o = -RsTM^ = ^)- 



(38) 



Equation (^6|) provides the connection between the mi- 
croscopic and selfconsistent coupling constants. Namely, 



1 



1 



self 



X' 



+ Rsv.ST{i^ = 0). 



(39) 



Hence, one can notice that the difference between the 
selfconsistent and symmetry restoring strength is related 
to the coupling of the field A to the kinetic energy opera- 
tor. Whereas the selfconsistent strength has to be calcu- 
lated using the prescription (||), the symmetry restoring 
one is defined as the inverse of the response function at 
zero frequency and guarantees that there must exist a 
spurious zero-energy RPA mode. 

One can easily notice that Rsv,5t{^ = 0) vanishes e.g. 
when we consider dipole distortions. The operators 
in this case has the form: 



Aw = VQio • p Pz 

Ml ^ V((5ii - Qi-i) - p^Px 

^fi = iV((3ii + Qi-i) - p^Py 
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